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1. Introduction 


Today, one common type of problem being solved on computers is the 
numerical integration of the initial value problem 

y'~ f(x,y) y(x Q ) = y Q . (l.l) 

Although there exist specialized techniques for solving such problems, 
the two most commonly used are the single step (e.g., Runge-Kutta) and 
multistep (e.g., Adams) methods. The major advantage of the multistep 
methods over the single step methods is that fewer functional evalua- 
tions are usually required per integration step. For this reason, the 
authors will limit the discussion that follows to predictor-corrector 
schemes „ 

For the individual who wants a solution to (l.l), there are two 
important considerations when selecting a method: speed and accuracy. 

For the researcher who is developing a predictor-corrector procedure, 
there are three important considerations: speed, truncation error and 

stability. The basic factor for determining speed is the number of 
functional evaluations required per step. Accuracy is determined by 
the stability error and truncation error of the method. The Adams 
methods have been shown to be about as good as any of the existing 
predictor-corrector formulas. The major problem with these formulas 
is that interval modification and interpolation are difficult. 

Nordsieck [l] has devised a method equivalent to one of the Adams 


1 



methods which accomplishes these tasks simply and inexpensively 
(Appendix C). Also, Nordsieck's method has a built-in automatic 
starting procedure, a feature not available with the Adams formulas, 

2 . The Methods 

The derivation of Nordsieck's method is based on Taylor's theorem. 
The authors refer the reader to Nordsieck’s paper for a detailed de- 
scription of the method. Nordsieck stores the current values of the 
higher derivatives of a polynomial approximating the solution. Adams 
stores successive values of a polynomial approximation for f(x,y) at 
several backward points. Nordsieck saves the equivalent polynomial 
information as Adams but in a more convenient form for interpolation 
and interval modification. 

The formulas discussed in detail by Nordsieck are equivalent to 
the 5th order Adams -Bashforth predictor and the 6th order Adams -Moulton 
corrector. Descloux [2]. has derived a Nordsieck-type formula which is 
equivalent to the 6th order Adams-Bashforth predictor and the 6th 
order Adams -Moulton corrector. The authors will study these particular 
formulas because of their desirable stability and truncation error 
properties. The following terminology will be adhered to: the original 

formulas derived by Nordsieck will be called the standard Nordsieck 
formulas; the formulas derived by Descloux will be called the modified 
Nordsieck formulas. These equations written in terms of backward 

points are given in (2.1) and (2.2). 


2 



Standard Nordsieck 


predictor: y(x+h) = y(x) + -~q (l901f(x) - 2774f(x-h) 

+ 26l6f (x-2h) - 127lff(x-3h) 

+ 251f(x-4h)) 

( 2 . 1 ) 

corrector: y(x+h) = y(x) + (475f(x+h) + l427f(x) 

- 798f(x-h) + 482f (x-2h) 

- 173f(x-3h) + 27f(x-4h)) 

Modified Nordsieck 

predictor: y(x+h) = y(x) + (4277f(x) - 7923f(x-h) 

+ 9982f(x-2h) - 7298f(x-3h) 

+ 2877f(x-4h) - 475f(x-5h)) 

( 2 . 2 ) 

corrector: same as standard Nordsieck corrector 

The corresponding equations derived by Nordsieck and Descloux are given 
in Appendix A. 

The reason for studying the modified Nordsieck formulas is that a 
step by step check on the local relative truncation error can be made. 
No estimation of the local relative truncation error is available with 
the standard Nordsieck formulas. 

When selecting a predictor-corrector scheme, one must decide the 
number of corrections (and evaluations) to be made after doing the 
initial prediction. This is of importance since this choice affects 
the stability and truncation error of the method. If one iterates to 
convergence then only properties of the corrector influence the result. 
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Because it is generally considered that functional evaluations are the 
most expensive part of the predictor-corrector procedure, the authors 
will limit the number of functional evaluations to two per step and 
test only the following procedures: PEC, PECE and PECEC. (P stands 

for predict, E stands for evaluate and C stands for correct.) 

3- Truncation Error 

For the single equation y ' = f(x,y), we use the notation: 
y(x) = true solution of the differential equation 
y ( x ) - Nordsieck solution for P(EC) 11 "’ X or PE(CE) li 

6 n (x) = y n U) - y( x ) 

f n (x) = f(x,y n (x)) 
y (x) = predicted value of y 
f P (x) as defined by Nordsieck 
h = the stepsize . 

The Nordsieck iteration for one step from x to x*h is then 
y (x+h) = y (x+h) (3.1) 

1 ir 

y n (x+h) = y (x+h) + h^^x+h). - f P (x+h)). (3-2) 

n^2 

The error for the standard Nordsieck is 

6 2 (x+h) = 5 § §o fr 7 y^(xH-h) - (g§8 "J h 7 y^ 6 ^ (x+h) (x+h,y(x+h)) ^ 

+ 0(h 8 ) 
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§ n (x+h) = h 7 y^(x+h) + 0(h 8 ) n^3 

and the error for the modified Nordsieck is 

6 n (x+h) = h 7 y ( 7 } ( x+h ) + 0(h 8 ) n^2 . 

In the standard Nordsieck case we know 

f (x+h) - f P (x+h) = h 5 y^ 6 ^(x+h) + 0(h 6 ) n^l 

y n (x+h) - y^(x+h) = h 6 y^ 6 ^(x+h) + 0(h 7 ) n^2 

y 3 (x+h) - y g ( x +h) = (g^-) 2 h 7 y^(x+h) (x+h,y(x+h)) 

+ 0 ( h 8 ) . 

f 2 (x+h) - f x (x+h) = 2§i" h6 y^ S ^( x+h ) ( x+h, y (x+h)) 

+ 0(h 7 ) , 


( 3 - 4 ) 


(3-5) 


(3-6) 

(3-7) 


(3-8) 


(3-9) 


(7) 

No satisfactory measure of y exists. However in the modified 
Nordsieck case we have the information 


f ( x+h ) - f P ( x+h ) = h 6 y^(x+h) + 0(h 7 ) n>l (3.10) 

y Q (x+h) - y p (x+h) = h 7 y'* 7 ' > (x+h)+0(h 8 ) n^2. (3*11) 
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Thus, in this case, an approximation for 6^ (x+h) is 


6 n~ iHo (y n (x+h) -y* +h » 1122 • 


(3.12) 


For a system of m equations 


y^' - f^(x,y (x), ..., y m (x)) , i = 1, ,.,m , we use the notation: 


y. (x) - true solution of the differential equation 

y. (x) = Nordsieck solution for P(EC) n ^ or PE(CEl n 
"■Lyi ’ ' ' ' 


6 . (x) = y. (x) - y. (x) 

i,n' J i,n l 

f. (x) = f . (x,y (x), y (x)) 

y^ p( x ) = Predicted value of. y^ 
f^(x) as defined by Nordsieck. 


The error for the standard Nordsieck is 


6 i, £ (x+h) - -wm h7 y i <7, <*+h) 


-21U J Y y ( 6 }( X+ h) 

288/ n L y j " ] 


(6), . , ,y i-n (^h),.., y (x + h)) 

I ( x+h ) 2 ^ 2 

dy. 


+ o(h a ) 


(3-13) 


i — 1, . . . , m 


6 l,n (x+h) = S$3o h t 1 (7) (x+h) + 0(b») n*3 

i = 1, . . . , m 


(3-14) 
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and the error for the modified Nordsieck is 


6, Jx+h) -^h 7 y, (7) (x+h) + 0(h 8 ) n£2 


i,n’ 


i - 1, . . ., m. 


(3-15) 


Equation (3.12) now becomes 


6 . 

1. 


n 


863 

19950 


(y, J x+h ) 

LI 


y. (x+h)) n>2 . 


(3-16) 


4. Stability 

The necessary theoretical background on stability can be found in 
the papers by Krogh [ 3 ], Crane and Klopfenstein [4] and Chase [ 5 ]. By 
examining the eigenvalues {1} of the Jacobian matrix of f with respect 
to y, the stability of a method can be investigated. It will be as- 
sumed that the Jacobian matrix is completely diagonal! zable and the 
eigenvalues are approximately constant over an interval of 5h (standard) 
or 6h (modified) . The authors have some evidence that the results in 
this section hold for Jacobian matrices which are not completely 
diagonalizable . Consider the set fs : s = hX] . The values of the 
elements of this set determine whether or not the method remains stable . 
Briefly, a method is absolutely stable if errors decrease in magnitude. 
A method is relatively stable if errors do not grow more rapidly than 
the solution. 

Formulas given by Hall [6] allow the computation of the character- 
istic equations 
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p( z, s ) = 0 (4.1) 

for the standard and modified Nordsieck PEC, PECE and PECE.C algorithms 
(Appendix B). With the appropriate conditions [5>P-38o] applied to 
these equations, the absolute and relative stability diagrams have been 
computed (figures 1-5)* Because of the symmetry involved, only the 
upper half of the s -plane is given in each diagram. The heavy solid 
curves bound the absolute stability regions. The portion of these 
curves which appears to coincide with the imaginary axis is just very 
close to it. The area to the right of the dashed curves is the region 
of relative stability. In the right half s -plane, this region does 
not have a well-defined boundary. Examination of the five diagrams 
yields the following two results: 

1) the diagrams for both the standard and modified Nordsieck 

formulas decrease in area in the following order: PECE, 

PECEC, PEC; 

2) the standard Nordsieck PEC, PECE and PECEC diagrams are 
larger than the modified Nordsieck PEC, PECE and PECEC 
diagrams respectively. 

Although diagrams for different order Adams formulas have not been 
computed, the same general result is expected. 

An example will show that absolute stability will not always 
insure an accurate solution. The stronger condition of relative stability 
will be required. The definition of relative stability presented by 
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Krogh [3, p.377] is as follows: 


Example . 


Defn. For |s| sufficiently small, one of the 

roots (of the characteristic equation (4.1) ) 

s 

r^, the principal root, is approximately e . 

Any other root we label r g to indicate that 
it is extraneous. If each of the r satisfy 
|r e |^|e S [, and if the r g with magnitude of 

g 

e are simple, then a method, is said to be 
relatively stable . 

Problem: y* = - y , y(o) = 1 

Solution: y(x) = e X , X = - 1 

Procedure : standard Nordsieck PEC with constant s 


Absolute Value of Relative Error 
s = -.07 s = -.08 


3 

•571-5 

.116-4 

6 

.486-5 

.272-4 

9 

•524-5 

.120-3 

12 

•533-5 

.251-3 

15 

•547-5 

.252-2 

18 

•561-5 

.656-2 

21 

•574-5 

.620-1 

24 

.584-5 

.159-0 
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Inside the relative stability region (s = -.07), the standard Nordsieck 
PEC calculates a good solution. However, outside the relative stabil- 
ity region but inside the absolute stability region (s = -.08) the 
solution blows up. The standard Nordsieck characteristic equation has 
been solved for s = -.07 and s - -.08. Below, e and the roots are 
listed. Complex conjugates are omitted. 


s 

-.07 

-.08 

s 

e 

. 9323938 

.9231164 

r i 

•93239^0 

.9231166 

r 2 

-.9135050 

-.9632736 

r 3 

-.305U55 + .140721 

.3120278 + .14246211 

r 4 

-.0811852 + .4828491 

.089245 + .4942325i 


Since |s! is sufficiently small, in both cases r , the principal root, 

approximates e S closely. When s = -.07, |r 2 |, | r 3 | and ! r^ | are all 
s 

less than e as Krogh has required for relative stability. However, 
when s = -.08, jr | is greater than e S . The root dominates the 
solution of the error difference equation [7,p.293l and thereby causes 
the poor results . 

s 

In the right half s-plane, as the real part of s gets larger, e 
"moves away" from the roots of the characteristic equation. When the 

g 

principal root no longer approximates e , then the definition of rela- 
tive stability can no longer be applied. It is then possible that the 
h is so large that the truncation error will cause a poor solution to 
be calculated. 
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Figure 1. Standard Nordsieck PEC 




Figure 2. Standard Nordsieck PECE 



Figure 3 - Standard Nordsieck PECEC; modified Nordsieck PECE 






Figure 4. Modified Norasieck PEC 







5. Computational Considerations 

Although the stability diagrams for the standard and modified 
Nordsieck formulas have been computed, these diagrams are difficult to 
use in many test problems. The difficulty lies in the fact that an 
eigenvalue problem must be solved at each integration step. Even if 
the partial derivatives are known functions or can be closely approxi- 
mated, the resultant eigenvalues may not be accurate. Matrix norms [ 8 ] 
are available which allow the computation of an upper bound for the 
moduli of the eigenvalues. However, there is no assurance that these 
approximations will be reliable. For these reasons, it is desirable 
to find some way to insure stability at each step without having to 
consider the values of s. Unfortunately, there is no known technique 
available to accomplish this for all initial value problems. 

Even if the eigenvalues are known, it is the nature of the true 
solution to the differential equation that determines which values 
of s are essential to the stability of the problem. Furthermore, the 
type of stability that is appropriate to consider is also determined 
by the behavior of the true solution. In many cases, relative stability 
is too strong a condition to impose to control the stability error. 


6 . Test Problems 

Because of the lack of a minimal comprehensive set of test problems 

for comparing the efficiency of our algorithms, the authors chose the 
following examples without claim to the completeness of the selections: 
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3) y 1 -2xy , y(o) - 1 , x 0 = 0 , x^ = 10 , 

X 2 

with solution y(x) = e and X = 2x ; 

4) y' = -xy(l + y 2 ) , y(o) = 1, x = o, x mov = 10, 

u mcix 

with solution y(x) = (2e x - l)~* and. 



The initial h for all problems was 1. Total and successful functional 
evaluations do not include evaluations done in the starting procedure. 
Since no accurate measure of the truncation error is available for the 
standard Nordsieck formulas, the authors defined 

6 = h(f (x+h) - f^(x+h)) for PEC, 

2 i 

- h(f (x+h) - f P (x+h)) for PECE, (6.l) 

6 = h ( f ( x+h ) - f ' ^ ( x+h ) ) for PECEC . 

3 2 
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was 


A local relative truncation error test ( I 6 n /y n (x + h) I < e) 
applied to test problems 1, 3 and 4. An absolute error test 
( ! 6 | < e) was applied to problem 2. The input values of e for 
our test programs were 5-0 x 10 4 , 2.5 x 10 4 , 10 4 , 7«5 x 10 5 , 
q 0 x io“5 % o o 3 y .5 x 10 -17 . No stability requirements were 
imposed in the test programs. The test problems were solved on an 
IBM 360/91 using double precision arithmetic. 

The results for e = 10" s and lo " 13 are given in Tables l-i+ . 

In an effort to make meaningful comparisons, graphs were drawn for 
each test problem plotting the final accumulated relative error 
(absolute error for problem 2) versus total functional evaluations 
for each value of e. These graphs are not presented here. 

7. Alternative Interval Modification 

An alternative method for determining the stepsize for the modified 
Nordsieck formulas was developed as follows: Assume that one step has 

been taken from x to x + h. Impose the relative error condition 


n 


y n (x+h) 


< e 


(7.l) 


Since 6^ - h 7 y^(x + h) + ... then an "optimum" interval may 

be found using the relation 


e 


y n ( x +h) 



(7.2) 
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or 


opt 


h 7 e 


y n (x+h) 


n 


Using the approximation (3-12) for 6^ yields 


y n (x+h)| 


. _ h 19930 e 

opt ~ ^ 863 | y n (x+h) - y^(x+h) 


i_ 

7 


( 7 . 4 ) 


For m-dimensional systems, let h ^ l 16 the minimum of the right 
hand side of ( 7 . 4 ) over the m equations. A similar analysis will 
hold for this section if one imposes an absolute error condition 
in ( 7 . 1 ). The appropriate changes are easily made in ( 7 . 4 ) and in the 
following strategy table. Numerical results appear in Tables 5 
and 6 . 
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Strategy Table 



t is a constant which satisfies 0 < t < 1 




8. Final Remarks 


Examination of the graphs discussed in section 6 revealed that, for 
high accuracy requirements (input e < 10 9 ), the standard EEC was more 
efficient than the standard PECE and PECEC procedures and the modified 
PEC was more efficient than any of the other procedures tested. For the 
modified Nordsieck formulas, this result held true for both techniques of 
interval modification. Neither the Nordsieck strategy nor the proposed 
strategy was clearly superior to the other. Further investigations are 
needed to determine when single evaluation predictor — corrector methods 
are more efficient than multi -evaluation methods. 

In conclusion, the authors note that subroutines employing Nordsieck's 
method have been satisfactorily used by members of the Goddard Space Flight 
Center for the past three years. The ease of performing interval modifi- 
cation and interpolation are the particular features of this method that 
were found to be most desirable by scientific personnel at this installation. 
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Table 1. Problem 1, 



starting 

h 

relative error 

of y at x 
J max 

std. PEC 

e = lO -6 2“ s 

.220-5 

" PECE 

ii 2~ 4 

.656-6 

" PECEC 

tt 2~ 3 

.111-6 

mod . PEC 

ii 2~e 

.888-6 

" PECE 

it 2~ 3 

.141-4 

" PECEC 

ii 2" 3 

•574r5 

std. PEC 

e = 10“ 13 2" 11 

.168-11 

" PECE 

It 2*"^ 

.177-11 

" PECEC 

tt 2***^ 

.181-11 

mod . PEC 

n 2” 11 

.296-10 

" PECE 

tt 2~7 

.245-10 

" . PECEC 

tt 2 _7 

.217-10 


no. of 
integration 
steps 

no. of 
successful 
evaluations 

no. of 
total 

evaluations 

no . of 
•evaluatio: 
to start 

1607 

1607 

1607 

121 

1600 

3200 

3200 

177 

1597 

3194 

3196 

145 

2376 

2376 

2410 

151 

800 

1600 

1600 

181 

800 

1600 

1600 

181 

25611 

25611 

25611 

201 

25603 

51206 

51206 

337 

25603 

51206 

51206 

337 

6425 

6425 

6425 

251 

6403 

12806 

12806 

341 

6403 

12806 

12806 

341 



Table 2. Problem 2. 




starting 

h 

absolute error 

of y, at x 

max 

no . of 
integration 
steps 

no. of 
successful 
evaluations 

no. of 
total 

evaluations 

nq. of 
evaluations 
to start 

std. PEC 

m 

II 

h-* 

O 

l 

ON 

2" 6 

.221-6 

327 

327 

327 

121 

" PECE 

fl 

2~ 4 

.653-7 

320 

640 

640 

177 

" PECEC 

11 

2 4 

.173-7 

320 

640 

640 

177 

mod . PEC 

?r 

2 6 

.534-6 

330 

330 

333 

151 

" PECE 

1! 

2 3 

.112-3 

83 

166 

166 

181 

" PECEC 

II 

2 3 

.860-4 

83 

166 

166 

181 

std. PEC 

in' 13 
e= 10 

2~ U 

.737-13 

5131 

5131 

5131 

201 

" PECE 

ii 

2~ 9 

.808-13 

5123 

10246 

10246 

337 

" PECEC 

ii 

o-9 

.848-13 

5123 

10246 

10246 

337 

mod . PEC 

H 

2~ U 

.343-11 

1305 

1305 

1305 

251 

" PECE 

ii 

2" 8 

.246-11 

1288 

2576 

2576 

381 

" PECEC 

n 

2~ 7 

.199-11 

1283 

2566 

2566 

341 



Table 3. Problem 3 


starting 

h 


std o 

PEC e = 

10 - 6 

2" 5 

ft 

PECE " 


2 -4 

n 

PECEC " 


2^ 

mod . 

PEC " 


2" 5 

f 1 

PECE " 


2 -3 

1! 

PECEC " 


2" 3 

std . 

PEC ® = 

10- 13 

2' 9 

If 

PECE 


2 8 

H 

PECE.G " 


2' 8 

mod . 

PEC " 


K) 

l 

00 

It 

PECE " 


2 7 

M 

PECEC " 


2‘ 7 


relative error 

of y at x 
J max 

no. of 
integration 
steps 

.615 -5 

1501 

.185 -5 

1501 

.178-6 

1499 

.470-3 

527 

.683-4 

518 

.113-3 

528 

.217-11 

22822 

.169-11 

22791 

.147-11 

22791 

.597-10 

5987 

.707-10 

6125 

.823-10 

6130 


no. of 
successful 
evaluations 

no . of 
total 

evaluations 

1501 

1506 

3002 

3010 

2998 

3008 

527 

543 

1036 

1044 

1056 

1066 

22822 

23015 

45582 

46052 

45582 

46052 

5987 

6117 

12250 

12256 

12260 

12266 


no. of 
evaluations 
to start 

105 

177 

177 

131 

181 

181 

169 

305 

305 

191 

341 

341 



Table 4. Problem 4 . 

starting 

h 

relative error 

of y at x 
J max 

no . of 
integration 
steps 

no. of 
successful 
evaluations 

no. of 
total 

evaluations 

no. of 
evaluation 
to start 

std. PEC e 

- IQ' 6 

2' 5 

.240-5 

745 

745 

773 

105 

" PECE 

II 

-4 

.162-5 

680 

1360 

1448 

177 

" PECEC 

II 

2 

.448-7 

709 

1418 

1418 

177 

mod . PEC 

It 

2~ 5 

.621-6 

1218 

1218 

1233 

131 

" PECE 

II 

2~ l 

.786-4 

318 

636 

685 

121 

" PECEC 

M 

t 

CNJ 

.148-4 

309 

618 

644 

241 

std. PEC e 

- 10- 13 

2' 9 

.114-12 

9942 

9942 

10481 

169 

" PECE 

It 

2 8 

.421-12 

10303 

20606 

20972 

305 

" PECEC 

It 

2~ 8 

.524-12 

10302 

20604 

20972 

305 

mod , PEC 

It 

2~ 8 

.528-10 

2717 

2717 

2805 

191 

" PECE 

II 

2' 1 

.322-10 

2863 

5726 

5790 

121 

" PECEC 

It 

2“ 7 

.277-10 

2807 

5586 

5614 

341 



Table 5 

. Modified 

-13 

Nordsieck with e = 10 , 

t = .95 



Problem 

Procedure 

* 

no. of 
integration 
steps 

no. of 
successful 
evaluations 

no. of 
total 

evaluations 

1 

PEG 

.465-9 

4185 

4185 

4185 

1 

PECE 

.352-9 

4166 

8332 

8332 

1 

PECEC 

.289-9 

4163 

8326 

8326 

2 

PEC 

.796-10 

817 

817 

817 

2 

PECE 

.507-10 

800 

1600 

1600 

2 

PECEC 

.369-10 

795 

1590 

1590 

3 

PEC 

.179-9 

4279 

4279 

4324 

3 

PECE 

.315 -9 

4270 

8540 

8630 

3 

PECEC 

.373-9 

4274 

8548 

8638 

4 

PEC 

.198-9 

1955 

1955 

2002 

4 

PECE 

.147-9 

1953 

3906 

4008 

4 

PECEC 

.121-9 

1951 

3902 

4004 


■*The values in this column correspond to the values in the second labeled columns of 
Tables 1-4, The starting procedure was unchanged for these test runs. 



Table 6. Modified Nordsieck with e = 10 


9 


t 


95 


-6 


Problem Procedure 

* 

no. of 
integration 
steps 

no. of 
successful 
evaluations 

no. of 
total 

evaluations 

1 

PEC 

o 232-5 

2434 

2434 

2592 

1 

PECE 

.695 -3 

460 

920 

922 

1 

PECEC 

.114-3 

453 

906 

908 

2 

PEC 

.626-5 

323 

323 

348 

2 

PECE 

.973-4 

82 

166 

166 

2 

PECEC 

.857-4 

84 

168 

168 

3 

PEC 

.150-2 

390 

390 

422 

3 

PECE 

.254-3 

377 

754 

826 

3 

PECEC 

.296-3 

382 

764 

836 

4 

PEC 

.847-5 

1244 

1244 

1337 

4 

PECE 

.297-3 

220 

440 

512 

4 

PECEC 

.331-4 

222 

444 

502 


*The values in this column correspond to the values in the second labeled columns of 
Tables 1-4 „ The starting procedure was unchanged for these test runs. 



Appendix A 


\ 


Standard Mordsieck Formula 

y(x+h) = y(x) + h(f(x) + a(x) + b(x) + c(x) + d(x) + Jgg (f(x+h) - f^Cx+h)))^ 
f p (x+h) = f(x) + 2a(x) + 3b(x) + 4c(x) + 

a(x+h) = a(x) + yoi-x.) + 6c(x) + 10d(x) + (f(x+h) - f P (x+h)), 

b(x+h) = b(x) + 4c(x) + 10d(x) + (f(x+h) - f P (x+h)), 

c(x+h) = c(x) + 5d(x) + (f(x+h) - f P (x+h))j 

d(x+h) = d(x) + (f (x+h) - f P (x+h) ) . 

Modified Nordsieck Formula 

y(x+h) = y(x) + h(f(x) + a(x) + b(x) + c(x) + d(x) + e(x) 

+ (f(x+h) - f P (x+b))), 

f P (x+h) = f(x) + 2a(x) + 3b(x) + 4c(x) + 5d(x) + 6e(x), 

a(x+h) = a(x) + Jb( x ) + 6c(x) + 10d(x) + 15 e(x) + (f(x+h) - f P (x+h)), 

b(x+h) = b(x) + 4c(x) + 10d(x) + 20e(x) + g (f(x+h) - f P (x+h)), 

c(x+h) = c(x) + 5<3-^ x ) + 15 e(x) + ^ (f(x+h) - f P (x+h)), 

d(x+h) = d(x) + 6e(x) + ^ (f(x+h) - f p (x-fh)),, 

e(x+h) = e(x) + ~q (f(x+h) - f P (x+h)). 
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Standard 


PEC: z 6 


PECE: z 5 


PECEC : 


Appendix B 

Characteristic Equations 


Nordsieck 


,, 4277 n r 2641 „ 

(1 + iss s ’ z= + w 5 z 


s z ! 

720 



s z 




= o 


- (1 + l5o 3 


36119 

+ 41472 ' 


z 4 + 


(133 s 
l 25o 


+ 


26333 s2 
20736 


M s 2) 
I726 ’ 




I2103 


“ ( i^ 3 + 4 i4t 2 s2) " 0 


(241 s 
l 720 


'4829 2 \ 
14-72 ' 


z 3 + 




69331 

4i472 


, 3 18221 _ 9025 2 

C lio 3 + 2755 H 3 5 z + W& 3 - 


0 
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Appendix B, Characteristic Equations, continued 


Modified Nordsieck 


PEC 




10 


+ s z 3 lg , 67 sz 2 + ^&sz--^s-0 

T 20 s z ‘ T5o s 2 + 258 s 2 8 u 


PECE: Same as standard Words ieck PECEC 


PECEC: z 7 - (1 + s + |2| s 2 ) z 6 


7 <S§ s + if sS > z5 - <# 5 + Ws sZ) 

* a -<& 


025 


/63175 2n 

a ) z - ^45 


= 0 
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Appendix C 


Modified Nordsieck Interval Modification 


Replaces h 
Replaces y 
Replaces f 
Replaces a 
Replaces b 
Replaces c 
Replaces d 
Replaces e 


Reversal 

Increase 

Decrease 

- h 

gh 

h/g 

y 

y 

y 

f 

f 

f 

- a 

ga 

a/p 

b 

p 2 b 

b/P 2 

- c 

p 3 c 

c/p 3 

d 

g 4 d 

d/p 4 

- e 

g 5 e 

e/p 5 


Modified Nordsieck Interpolation 


Suppose z is between x and x + h 


Let of = 


z-x 

h 


then y(z) = y(x) + h(of (x,y(x) ) + cPaix) + or^fx) + c/Mx) 


+ cr'd(x) + o^eix)) 


f(z,y(z)) = f (x,y(x) ) + 2ct?a(x) + Jc^bfx) + ka 3 c(x) 


+ 5 a 4 d(x) + 6oPe (x) . 


For standard Nordsieck case, omit all e(x) terms in the above formulas. 
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